Bulletin of Electrical Engineering and Informatics 
Vol. 13, No. 1, February 2024, pp. 233~246 
ISSN: 2302-9285, DOI: 10.11591/eei.v1311.5371 o 233 


BER estimation for STBC-MC-DS-CDMA-4 antennas system 
by varied wavelet-carriers features via AWGN-flat channels 


Nader Abdullah Khadim!, Ali Jawad Kadhim Jawad Alrubaie?, Osama Qasim Jumah Al-Thahab* 
'Medical Instrumentation Techniques Engineering, Al-Mustaqbal University College, Hillah, Iraq 
College of engineering, University of Babylon, Hillah, Iraq 
Department of Electrical, College of Engineering, University of Babylon, Hillah, Iraq 


Article Info 


ABSTRACT 


Article history: 


Received Nov 29, 2022 
Revised Apr 17, 2023 
Accepted May 10, 2023 


This paper improves the bit error rate (BER) of the modern communication 
system by taking into account the effect of the wavelet shape and the number 
of carriers on the performance of the space time block code multi-carrier 
direct sequence code division multiple access (STBC-MC-DS-CDMA). Here 
the transmitter is moved with speeds 2 km/hr, 45 km/hr and 100 km/hr via 


Rayleigh flat fading channel. Here, 2 antennas are employed at the receiver 


to mitigate the multipath signal influence. The system’s orthogonal 
frequencies are generated using Haar, Daubechies 4, Symlets 4, Cohen- 
Daubechies-Feauveau 1.1 with 9.7 and B-spline 3. The number of used 
carriers is 128, 512, and 1,024. Quadrature phase shift key (QPSK) is used 
with cyclic prefix 1/16 and a bandwidth of 20 MHz. traditional fast fourier 
transform (FFT) system is compared to the proposed discrete wavelet packet 
transform (DWPT) to show the BER enhancements. The space-time block 
coding (STBC) is used to enhance system cabability in error corection. The 
proposed system shows significant improvement in BER, so that, it reaches 
to the same BER when using FFT but with less signal to noise ratio (SNR), 
which interns reduce the power consumed within the system and the cost. 
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1. INTRODUCTION 

The recent advancements in data network gave rise to many applications that are involved with 
almost every aspect of life ranging from top secret military operation to merely sending greetings for a 
birthday [1], [2]. These networks are becoming the backbone of societies because people are relying on them 
to spread news all over the world [3]-[5]. Even the elections consider social networks as a reliable source of 
information to study the citizen’s point of views regarding a certain party or a candidate. Presently, these 
networks are growing to be self-sustained economic systems and a pillar of modern societies [6], [7]. 

To satisfy the higher bandwidths, the engineers could follow the traditional way by providing the 
extra bandwidths for the networks. Although this might be satisfactory it suffers from by unavoidable 
limitations. First limitation is the scares affordable free bandwidth because the radio frequency (RF) 
spectrum is highly congested. The second limitation is the bandwidth over a wide range will increase the 
influence of nonlinearity and deformations effect on the transmitted signals [8]-[11]. Differnt approaches 
will invented to modulate the data that used narrower bandwidth with higher data rates, like the multi-carrier 
transmission (MC) or its modern version the orthogonal frequency division multiplexing (OFDM) [12], [13]. 
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Khan et al. [14] suggested a system of using discrete wavelet packet transform (DWPT) in OFDM 
instead of discret fourier transform (DFT), then compared the performance between them based on bit error 
rate (BER) values and their corresponding signal to noise ratio (SNR), while Soni et al. [15], used a discret 
wavelet transform (DWT) instead of fast fourier transform (FFT) in the OFDM system and test the proposed 
system by sending various images through different channes. They found that, the DWT is more reliable and 
has less SNR as compared to FFT. The use of these techniques is very important in many application as 
Roy et al. [16] depictes in his paper, so that a comparison is made between the FFT, discrete cosine transform 
(DCT), DWT and wals-hadarmard transform (WHT) in many applications like electrocardiogram and 
photo-plethysmography signals. This paper aims to improve the BER of the code division multiple access 
(CDMA) multicarrier system by adding a DWPT instead of FFT as a suggestion to the system. Here the 
effected thing is reducing the SNR of the proposed system so that it reaches to the same BER but with low 
power, which consequently reduce the complexity and cost. 


2. MC AND OFDM SYSTEMS 

The idea of OFDM is based on dividing high bit rate stream of data into low bit stream segments. 
These segments are then modulated using orthogonal subcarriers to ensure neither overlapping between 
adjacent. These sub-streams are then combined together and transmitted using a main carrier before 
transmission. The block diagram of typical OFDM system is shown in Figure 1. Due to the fact that the 
resultant subs-terms have low rate with a low bandwidth requirement, then the OFDM requires less channel 
compared with the classical system. The fact that the subcarriers are orthogonal also allows the channels to 
overlap without interference because the orthogonality ensures complete separation among channels [12]. 
The carriers and pilots illustration are shown in Figure 2, so that, Figure 2(a) depicts the user and pilot 
subcarriers and how the packets are performed, while Figure 2(b) states the effect of orthogonality on the 
subcarrier [13]. 


Serialto | _ cP 
Parallel ‘Removal 


Figure 1. OFDM block diagram 


Consider a multi-carrier (MC) OFDM system with N subcarriers and each subcarrier fẹ is used to 
modulate są sub-frame, then s(t) transmitted signal is (1): 


S(t) = Liter Skexp G2mk fit) (1) 
Here fp is equally spaced subcarriers and f, = BW /N, BW is the bandwidth. 


Let T to be 2/BW, then according to Nyquist theorem each sample occurs at T/2 or 1/BW. 
Therefore, any sub-stream subcarrier will exist at time mT then (1) becomes: 


s(mT) = YN_, s,exp (j2 =) (2) 
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In (2) is exactly the inverse discrete fourier transform (IDFT) or inverse fast fourier transform 
(IFFT). At the receiver, an FFT circuit can be used to retrieve the transmitted signals from the incoming 
samples. This is quite clear from Figure 1, where the symbols are modulated, possibly using quadrature 
amplitude moduation (QAM) [13]. 
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Figure 2. OFDM idea; (a) the subcarriers assignment to users with pilot carriers and (b) OFDM signal 
spectrum 


3. DISCRETE WAVELET PACKET TRANSFORM DWPT 
In (2) can be viewed as a running filter with coefficients exp (j 2T =) or in other words a low pass 


filter bank. The frequency location for each subcarrier fẹ occurs at every 2* x fọ, where fo is the fundamental 
frequency and fẹ are its harmonics. This means that the power spectral density (PSD) will be widely spread 
over the channel as the sinusoidal functions extends from -œo to +00. This consequently reduces the system 
spectral efficiency. In order to solve this problem, the analyzing function exp (j 21 =) can be replaced by a 


more compact or spectral efficient function to ensure higher PSD concentration [14]. The DWT can ensure the 
achievement of this goal. Hence; by replacing DFT or FFT with DWT, the spectral efficiency can be increased 
due to time limited period of DWT functions resulting in a more reliable and higher signal quality [15], [16]. 
The wavelets are orthonormal functions such that [17], [18]: 


1 l=k 


0 l#k 6) 


(PE, P) = S POROA = { 


Where Y(t) is the mother wavelet, and (-) is inner product. The p vanishing moments can be found as (4): 


co 


J, t*¥(Q)dt = 0 forO <k <p (4) 
The vanishing moments describe the speed of decay of the wavelet function. So, the DWT output is: 


Xw(s,r) = f°. x(t)¥,,(t)dt (5) 
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Where s,r is the scale and translation. This means that the wavelet can capture the details of x(t) in both 
time and frequency unlike FFT, which focuses only on the frequency domain. The translation means the 
wavelet is shifted in time while the scale is related to the width of the finite wavelet function. The scale and 
translation are shown in Figure 3 [19]. Assume an S samples discrete signal x[n] (n = 0,1,-:-,S — 1) want to 
apply first level DWT to it, then: 


yin] = Ek=-o x[k]Go[2n — k] + Veo x[k]G [2n — k] = Yow [N] + Ynign [n] (6) 
Scaling Translation 
W(t) V(t) 
a= 1 
t t 
(21) Y(t-t) 
a= 1/2 


Figure 3. The scale and translation of wavelet 


In (6), the Yiow is the signal coming from the low pass filter (LPF) bank with its impulse response 
Go[2n— k], while, Ynign is the signal coming from the high pass filter (HPF) bank with its impulse response 
G,[2n-k]. It is noted that the out coming signal has double the sampling rate n as the input signal. Thus, a 
decimation circuit by 2 (2) is required to switch back to the original sampling rate. Figure 4 shows the filter 


bank of the multi-level DWT [20]. The DWT functions Gp and G, suggest that the incoming signal will be 
x 
levels. The reconstruction functions Hy and H, are the inverse for the analyzing functions Gp and G4. In 
reconstruction filter, the sampling rate is up by a factor of 2 (2) to reverse the decimation process. The DWT 
can be modified to be identical to FFT and further increases the final signal PSD concentration. The 
modification comes from decomposing Yhign in the same manner y} at every stage. This analysis is called 


DWPT. The filters used for DWPT is shown in Figure 5 [21], [22]. 


transformed to orthogonal frequencies in the rate of = where k = 1, 2,:--,K which is the number of analysis 
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Figure 4. Filter banks and LPF response 
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Figure 5. DWPT analysis filter 


To generate the orthonormal subcarriers by using inverse discrete wavelet packet transform 
(IDWPT), consider (6) and by using mallat’s pyramid algorithm (MPA) for multi-resolution analysis (MRA), 
in (6) can be written as (7), so that g(k) is called the wavelet coefficient function and is related to h(.) (10). 
Then for a scaling function of order N, in (11) can be depicted. In addition, the dilation and wavelet 
equations can be written as seen in (12) and (13): 


W(t) = Eeo Ay (KIO), (E K) + LiL jy Eke- Bj (KY (2™t — k) (7) 
aj(k) = Em h(m — 2k)aj41(m) (8) 
bj(k) = Xm g(m — 2k)bj+ı (m) (9) 
g(k) = (-1)*h(1 =k) (10) 
g(k) = (D*AN -1 - k) (11) 
DCE) = V2 Zkzo A(k)@(2t — k) (12) 
Wt) = Yat (-1)*h(k)@(2t +k —N +1) (13) 


Therefore, by using MPA and extending it to the DWPT, series of x[n] vector can be analyzed using the 
following transformation matrix for L number of stages [23], [24]. 


hO) h(1) A(Z) h(3) ©  h(L-1) 0 0 i 0 0 

0 0 h(0) hA) A) B) © h(L-1) 0 T 0 

h(4) si h(L—1) 0 0 a 0 h(0) h() h(2) h(3) 
_|A(2) G3) as h(L-1) 0 0 a 0 0 RO) ACI) (14) 
Jg I g2) g(3) ve g(b=1) 0 0 m 0 0 

0 0 sO sO 9@) I® v -D O = o 

g4) ei g(0) 0 i 0 0 g(0) g(1) g(2) 9(3) 

g(2) g(3) whe g(L-1) 0 0 s 0 o g@) g0) 
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By substituting (11) in (14), the transformation matrix becomes: 


h(0) h(1) h(2) h(3) 5 h(L-1) 0 0 "i 0 0 
0 0 h(0) h(1) h(2) h(3) ~ AL=-1) 0 n 0 
h(a) ` A(L—1) 0 0 ` 0 h(0) h) h(2) hG) 
au| t h(3) = h(L—1) 0 0 N 0 0 h(0) h1) (15) 
“lh(L-1) -A(L—2) h(L-3) 2 hA)  —-h(0) 0 0 2 0 0 
0 0 h(L—1) -h(L—2) h(L-3) i h(1) —=h(0) 0 A 0 
h(l- 5) i h(0) 0 S 0 0 R-i) =A= n=) A= A 
h(L = 3) x, h(1) —h(0) 0 0 $ 0 0 h(L—1) -h(L-2) 


Using the orthonormal property, the reconstruction matrix is simply the transpose of the analysis function. 
Hence, at the construction side the following matrix is used (16): 


h(0) 0 0 0 0 h(2) A(L—1) 0 0 0 h(1) 
h(1) 0 0 0 0 : : 0 0 0 —h(0) 
h(0) 0 0 0 AL- AQ) hl-1) 0 0 0 
h(L—-1) R) 0 0 0 0 —h(0) 0 0 0 
0 : h(0) : 0 0 0 Ad) AO : 0 
0 h(L-1) RA) 0 0 0 0 —h(0) : 0 0 
| 0 0 : 0 : 0 0 0 —h(1) 0 : 
Ph ho 0 h(L - 1) 0 0 0 0 0 h(0) 0 0 (16) 
0 0 0 0 0 0 0 0 0 0 
: 0 0 0 0 0 : 0 0 0 
0 : 0 h(0) 0 : 0 : 0 AL- 0 
0 0 h1) 0 0 0 0 : hA(L - 1) 
0 0 0 ve : h(0) 0 0 0 0 —h(1) 
0 0 0 “ A(L-1) AC) 0 0 0 0 h(0) h(2) 


4. OFDM AND CDMA SYSTEMS 

The OFDM can ensure a complete separation for overlapped channels due to the orthogonality 
property. This property is tightly adhered to the location of subcarriers. If the location of these subcarriers are 
dispositioned, then the orthogonality will fall apart and the system will experience the inter channel 
interference (ICI) [24]. The ICI will reduce the system performance due to the decrease in the signal quality. 
Therefore, certain modifications should be taken, and one of them is using CDMA [25]. The CDMA is built 
upon the spread spectrum system (SSS) technology by which it allows all the users to use the channel all the 
time. A typical CDMA is shown in Figure 6 [26]. In the CDMA, each user is assigned by a pseudo 
noise-sequence (PN-sequence). The word pseudo comes from the fact that the sequence will repeat itself after 
a period of time depending on the linear feedback shift register (LFSR) constraint length. The generated PN 
pulses are almost random. The auto correlation function (ACF) of the PN sequence has the property given 
in (17) [27]. 


T=0 


T#0 (17) 


ACF(t) = | PN(t)PN(t + t)dt = 2 1 


Where N is the maximal length and equals to 2”-1 and n is the LFSR constraint length. 

When the PN sequence is aligned with itself then the ACF will have 1 value otherwise the sequence 
will cancel the incoming data. Unfortunately, the PN-sequence has poor cross correlation function (CCF), so 
gold code is used. If the sequences are chosen carefully, then the CCF will have only three states, which are: 


n+1 
{qt 2 200°} eq) = f +22 fornodd (18) 
, , = n+2 
N ye 1+22 forn even 


At the receiver side, the incoming streams are multiplied by the same gold code, then passed through a 
coherent detector as seen in Figure 7 [26]-[28]. Let the transmitted OFDM signal to be s(t), then the 
received signal y(t) is seen in (19): 


y(t) = s(t)e(t) + NÇ) (19) 


where C(t) is the spreading code at the transmitter side, N(t) is the impairments from the channel like noise 
or jamming or interference or multipath. 
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NCE) =n(t) +J (t) + Xita s(t — et — 7) (20) 
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Figure 6. A typical CDMA system block diagram 
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Figure 7. A detailed spread spectrum demodulator 


Where n(t) is the additive wide gausian noise (AWGN), J(t) is the jamming signal, t is the time delay due 
to multipath and M is the number of multipath signals. After multiplying y(t) by c(t) then output signal b(t) 
of the coherent detector for the input s(t) is: 


b(t) = = |" s(t)é(t)e(t)dt + f.” N(t)e(t)dt] (21) 


5. FREQUENCY FLAT FADING CHANNEL AND SPACE TIME BLOCK CODING 
The practical channels are usually imperfect and thus are described using its impulse response. assume 
a frequency selective fading channel is h of length L then the received signal r(t) is given in (22) [29]: 


r(n) = EE} h(k)s(n— k) (22) 


Where s(n) is the transmitted signal, in (22) shows that the channel with such channel impulse response (CIR) 
which will cause inter-symbols interference (ISI) or in the case of OFDM will be worse which is inter-block 
interference (IBI) or ICI. One of the methods that used to overcome the fading effect is the multiple-input and 
multiple-output (MIMO)-space time block coding (STBC) [30]. MIMO is a technique by which the symbol is 
transmitted through different paths in order to null the CIR. Here Alamouti’s STBC is used to avoid the 
continuous need to update the CIR [31]. 


BER estimation for STBC-MC-DS-CDMA-4 antennas system by ... (Nader Abdullah Khadam) 


240 o ISSN: 2302-9285 


The simplest form of Alamouti scheme is the multiple-input single-output (MISO) of 1x2 which 
means that the number of receiving antennas is 1 and the transmitting antennas is 2 as in Figure 8. Let two 
symbols x, and x, to be transmitted by this MISO over two channels of CIR or sometimes called channel 
state information (CSI) of h4 and hz, then the received symbols are: 


nO ha) 4] +r) (23) 


Tx antenna 1 l Rx antenna 


Path 1 >» 


path 252 
Tx antenna 2 


Mobile 
phone 


Figure 8. A 2x1 single input multiple output (SIMO) system 


At the second instance, the first antenna will transmit x3(t) and the second antenna will transmit —x{(t), so 
the received symbols at the second instance are: 


yo(t) = (hr ha %2,] +o (24) 
Taking the conjugate of (24) and arranging the result to look like: 
* * * x. * 
yi = (ny ife] +5 25) 


Where n,(t) and n2(t) are the AWGN in the first and second transmission instances respectively. 
Combining (24) with (25), then it can be shown that: 


Ole allo lee eB 
tet Y= oha = [is] co = [pi] ana = [Rae] en ¥ = Ee cal [GS] + 


Here, the Alamouti orthognal space time block code (OSTBC) can be extended from MISO to MIMO easily 
[32], [33]. Hı and H; is depicted in (27) and so on for H4. 


H, = a d and H; = 


Hy | 
Thy hi 


-H; Hi Gy 


6. RESULTS AND DISCUSSION 

The suggested system combines all the above concepts and integrate them into a single 
communication system as shown in Figure 9. First, the signals are mapped using quadrature phase shift key 
(QPSK) scheme. An OSTBC is added then the resultant stream is fed to IDWPT section to generate the 
orthogonal subcarriers. These pilots are used for synchronization and estimating the CSI of the channel in 
order to provide the necessary equalization if needed. This s erial stream is then fed to a MIMO of 2x2 
antennas then send over the channel. 

The transmitter is assumed to be a mobile travelling at speeds 2 km/hr, 45 km/hr, and 100 km/hr 
creating maximum doppler shift (MDS) of 10.7 Hz, 241.7 Hz, and 537 Hz respectively. The channel is flat 
fading rayleigh channel with AWGN. The number of transmitted bits per test is (102,400 bits) coming from 
(100 packet*64 characters *16 bit/char). The OFDM physical layer settings are shown in Table 1. 
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Table 1. OFDM physical layer simulation settings 


Parameter Mobile MC-DS-CDMA Scalable OFDMA-PHY 
FFT or DWT size 128 512 1024 2048 
Number of used data subcarriers 64 180 360 720 
Modulation types QPSK 
Cyclic prefix 1/16 
Channel bandwidth (MHz) 20 20 20 20 


Here, the mother wavelets utilized are as follows: Haar, Daubechies=db4, Symlets=sym4, 
Cohen-Daubechies-Feauveau=cdf, and B-spline 3=bs3. The results are as follows: 

a. MDS 10.7 Hz, the BER are shown in Figure 10 (see in Appendix) 
b. MDS 241.7 Hz, the BER are shown in Figure 11(see in Appendix) 
c. MDS 537 Hz, the BER are shown in Figure 12 (see in Appendix) 

From Figure 10, it can be seen that to obtain a 10° BER, the SNR must be at least 28 db, while by 
using the proposed system it decreaded to about 5 dB when using 512 subcarriers, and so on for the athors. 
The percentage improvent is about 80 percent, which is very good. The big difference between the results of 
the proposed system and the traditional method beyond to the excellent orthogonality of DWPT over FFT, 
which reduces the effect of ICI and ISI on the transmited data. 

A comparison is made when motor driven systems (MDS) equal to 10.7 Hz [30], they reach to 
10° BER (128 subcarrier) at 8 dB SNR in AWGN, while the proposed system reaches the same value at 5 db. 
In addition, Ali et al. [34] for the same sub carrier reches to 7 dB SNR, while at 1,024 subcarriers the SNR is 
8 dB and the proposed system of this work gives 7 dB. In the same context, Kusumawardhania et al. [35] 
obtain SNR of about 7.5 dB over AWGN channel when using a multi-level coding CDMA. This comparison 
gives the pereference of the proposed system over different algorithms used to enhance the BER evaluation. 


7. CONCLUSION 

Here, IDWPT/DWPT has an excellent performance as compared to the IFFT/FFT in all tests. This 
justifies the assumption that the DWPT has higher concentrated signal PSD. But in all cases, the BER 
increases with the increase of bit rate and this event is inevitable. The doppler shift has high impact on 
system performance as it increases due to the motion of the mobile station. The doppler shift changes the 
locations of the orthogonal subcarriers which derifts the orthogonality that depends on the difference between 
adjacent subcarriers, although, the pilots were transmitted to resynchronize the system but even at certain bit 
rate and motion speed they lose their effectiveness. Neither the STBC nor the CDMA are now able to track 
and fight of the doppler shift effect. It is highly noticed that the wavelets have no close relation to the bit rate 
and inspite of that, there are some kinds of wavelets that gives the worst performance upon FFT like bs3. 
Therefore, the wavelet selection should be carefully done. 
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APPENDIX 
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Figure 10. BER performance of STBC DWT-mobile MC-DS-CDMA in AWGN-flat fading channel with 
128, 512, and 1,024 subcarriers respectively 
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Figure 11. BER performance of STBC DWT-mobile MC-DS-CDMA in AWGN -flat channel with 128, 512, 
and 1,024 subcarriers respectively 
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Figure 12. BER performance of STBC DWT-mobile MC-DS-CDMA in AWGN-flat channel with 128, 512, 
and 1,024 subcarriers respectively 
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